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Abstract 

We consider one-dimensional, locally finite interacting particle systems with two con- 
servation laws. The models have a family of stationary measures with product structure 
and we assume the existence of a uniform bound on the inverse of the spectral gap which 
is quadratic in the size of the system. Under Eulerian scaling the hydrodynamic limit for 
the macroscopic density profiles leads to a two-component system of conservation laws. The 
resulting pdc is hyperbolic inside the physical domain of the macroscopic densities, with 
possible loss of hyperbolicity at the boundary. 

We investigate the propagation of small perturbations around a hyperbolic equilibrium 
point. We prove that the perturbations essentially evolve according to two decoupled Burgers 
equations. The scaling is not Eulerian: if the lattice constant is n^^, the perturbations are 
of order then time is speeded up by n^~^^ . Our derivation holds for < /9 < i. The 
proof relies on Yau's relative entropy method, thus it applies only in the regime of smooth 
solutions. 

This result is an extension of ^01 and where the analogue result was proved for 
systems with one conservation law. It also complements where it was shown that 
perturbations around a non-hyperbolic boundary equilibrium point are driven by a universal 
two-by-two system of conservation laws. 

1 Introduction 

There are several results dealing with the perturbation analysis of hydrodynamic limits for in- 
teracting particle systems. In the landmark paper |1] the authors prove that for the asymmetric 
simple exclusion, in dimensions higher than 2, perturbations of order n^^ of a constant profile 
evolve according to a certain parabolic equation under diffusive scaling (time rescaled by n^, 
space by n). It is well-known, that under Eulerian scaling (time rescaled by n, space by n) the 
hydrodynamic limit leads to a hyperbolic conservation law (the Burgers equation), the pertur- 
bation limit gives the same equation with the Navier-Stokes correction. (For a survey on the 
microscopic interpretations of the Navier-Stokes equations see the end of Chapter 7 of [HI.) 

Motivated by T. Seppalainen investigated a similar problem in one dimension for the 
so-called totally asymmetric stick process. In JOl he proves that an 0{n~^) perturbation of 
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the constant profile is governed by the Burgers equation (even after the appearance of shocks) 
if time is rescaled by and space by n, where /3 G (0, ^) is a fixed constant. Independently, 
in the authors partially extend this result by proving that one gets universally the Burgers 
equation in the hydrodynamic limit for similar perturbations of equilibrium for a wide class 
of one-dimensional interacting particle systems with one conservation law. The models are 
not reversible and not necessarily attractive. The proof relies on H. T. Yau's relative entropy 
method, it only applies in the smooth regime of solutions and it only works for /? G (0, -g). It is 
conjectured that the result should hold for all (3 G (0, ^) even without the smoothness condition 
as in the result of [Tn] . 

This universal result may be explained by the following arguments. Under Eulerian scaling 
these systems admit in the hydrodynamic limit a hyperbolic conservation law of the form 

dtu + d^j{u) = o. (1) 

Taking a point uq with J"[uq) ^ simple (although formal) calculations yield that solutions 
of (^, with initial conditions which are small perturbations uq, are governed by the Burgers 
equation. See for the 'more precise' formulation. 

In the present paper we give an extension of the results of lllj for systems with 2 
conserved quantities. In J2| a general one-dimensional family of lattice-models was introduced. 
The models are locally finite interacting particle systems with two conservation laws which 
possess a family of stationary measures with product structure. In that paper it is shown (in 
the regime of smooth solutions) that in Eulerian scaling we get a hydrodynamic limit of the 
form 

dtu + d^^u,v) = 0, 
dtv + d^^{u,v) = 0, 

where (u, v) V and P is a convex compact polygon, the the physical domain (see © for the 
definition). We also note, that |H] gives the first major result about the Eulerian hydrodynamic 
limit of multi-component hyperbolic systems. In it was also shown that an Onsager-type 
symmetry relation holds for the macroscopic flux functions (see Lemma One of the 
consequences of this relation is that inside the physical domain D the pde is (weakly) hyperbolic, 
i.e. the Jacobian can be diagonalized in the real sense. Experience shows, that the limiting 
pde is strongly hyperbolic (the Jacobian has two distinct real eigenvalues) in the whole physical 
domain except some special points on the boundary dV. 

We consider perturbations of order around a constant equilibrium point {uo,Vo) G P, 
which is strictly hyperbolic. We prove that rescaling time by n^"*"^ and space by n the evolution 
of the perturbations are governed by two decoupled equations. (These are 'usually' Burgers 
equations, see the remark at the end of subsection 13. 11 ) This result agrees with the formal 
perturbation of the pde ((2) e.g. with the method of weakly nonlinear geometric optics (see 

131 El). 
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The reason for the decouphng of the resulting pde system is the strict hyperbohcity, basically, 
the two different eigenvalues (sound speeds) cause the equations to separate. In the paper J2] 
perturbation around a special non-hyperbolic point was considered in a similar setting, it was 
proved that in that case in the limit the evolution obeys a two-by-two system of conservation 
laws which cannot be decoupled. The treatment of that problem needs more complex tools 
than our proofs, sophisticated pde methods are used besides Yau's method. 

Our proof follows the relative entropy method using similar steps as (thus it only applies 
in the regime of smooth solutions), but it also heavily relies on the Onsager-type symmetry 
relation proved in ^2]. We assume the existence of a uniform bound on the inverse of the 
spectral gap, quadratic in system size, to be able to prove the so-called one block estimate. We 
do not deal with the proof of the spectral gap bound, but we remark that with the techniques 
of [3 one can get the desired gap estimates for a large class of systems. Our result holds for 
P € (0, ^). Assuming the stronger (but harder to prove) logarithmic-Sobolev bound we could 
get the result for /3 € (0, |). 

2 Microscopic models 

We consider the family of microscopic models investigated in ^21 ■ We go over the definitions 
and the important properties, for the details we refer the reader to the original paper. There 
are several concrete examples introduced in J2j, we do not list them here. 

2.1 State space, conserved quantities, generator 

Throughout this paper we denote by T" the discrete tori Z/nZ, n G N, and by T the continuous 
torus M/Z. We will denote the local spin state by we only consider the case when Q is finite. 
The state space of the interacting particle system is 

Configurations will be denoted 

w := {LOj)j£jrL G ri". 
The two conserved quantities are denoted by 

we also use the notations (j = C{uJj), r]j = rj{ijjj). We assume that the conserved quantities are 
different and non-trivial, i.e. the functions C,, rj and the constant function 1 on $7 are linearly 
independent. 

We consider the rate function riOxOxfixO— > . The dynamics of the system consists 
of elementary jumps effecting nearest neighbor spins, {ujj,ujj^i) — > (ujpUJj^-^), performed with 
rate r{ujj,LL!j+i;ujj,uj'j_^-^^). 

We require that the rate function r satisfy the following conditions: 
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(A) If r(c(;i, C(;2; 1^1, ^2) > then 

C(a;i)+CM = CK)+CK), .3^ 
r]{uJi) + r]iuJ2) = r]{uj[) +r]{u;'2). 

This means that C and ^ are indeed conserved quantities. 

(B) For every Z G [n min ^, n max ("] fl Z, € [n min r/, n max 77] fl Z the set 

is an irreducible component of fi", i.e. if tj' G ^^at then there exists a series of ele- 
mentary jumps with positive rates transforming io into u/. This ensures that there are no 
hidden conservation laws. 

(C) There exists a probability measure vr on which puts positive mass on each element of 
and for any wi, u)2, uj[, UJ2 £ 0, 

Q{uJi,UJ2) + Q(W2, W3) + Q{uJ3,Ll!i) = 0, 

where 

Q{uJi,uJ2) := > <^ -7 — —, — -r{uJi,uj2;uJi,uJ2) -r{uJi,uJ2;uJi,uj2)> . 
This condition will imply that the measure HjeT" ^ stationary for our process on $7". 

U)' U)" 

For a precise formulation of the infinitesimal generator on we first define the map Qj ' : 

_^ fQj. g^gj.y ^f^^lf (ZQJ (Z jn. 

uj' if i = j 

uj" if i = j 

if i / j,j + 1. 

The infinitesimal generator of the process defined on 17" is 

We denote by X^^ the Markov process on the state space O" with infinitesimal generator U' 
2.2 Stationary measures 

For every ^, r € M let G{9, r) be the moment generating function defined below: 

G{e,T) := log J2 e''^^'^^+"''^"^vr(tj). 



ej'" ^J.= <; a;" if i=i + l 
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We define the probability measures 



T^eA^) ■= ^(^) exp(0C(c^) + rr?(c^) - G(0, r)) (4) 

on 0. Using condition (Q, by very similar considerations as in P, |^ or one can show 
that for any 0, r € M the product measure 

is stationary for the Markov process Af" on 17" with infinitesimal generator L". We will refer to 
these measures as the canonical measures. Since Qj and ?7j are conserved, the canonical 
measures on il" are not ergodic. The conditioned measures defined on ^ by: 



are also stationary and due to condition satisfied by the rate functions they are also ergodic. 
We shall call these measures the microcanonical measures of our system. It is easy to see that 
the measure tt^tv does not depend on the values of 9,t. 

2.3 Expectations, fluxes 

Expectation, variance, covariance with respect to the measures ttq^ will be denoted by Eg 
Var0,^(.), Co^rgA■)■ 

We compute the expectations of the conserved quantities with respect to the canonical 
measures, as functions of the parameters 9 and r: 



u{e, r) := EeAC) = Yl C(^)^sA^) = deG{e, r) = Gq, 
v{e,T) := Ee,r{v) = Y.v{^)^eA^) = drG{e,T) = Gr. 



We will usually note partial derivatives by using the respective subscripts, as long as it does 
not cause confusion. Elementary calculations show, that the matrix-valued function 

ug ^r\( Gee Ger \ =,G{9,r) 

Vg Vr J V ^Sr Grr J 

is equal to the covariance matrix Covg^ridv) ^^^d ^ consequence the function {9,t) i— 
{u{9,t),v{6,t)) is invertible. We denote the inverse function by {u,v) {9{u,v),t{u,v)). 
Denote by {u, v) S{u, v) the convex conjugate (Legendre transform) of the strictly convex 
function (6*, r) i-^ G{9,t): 

S{u, v) := sup (u9 + VT- G{9, r)) , (5) 
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and 

V := {{u,v) eR+ xR: S{u,v) < 00} (6) 

= co{(C(w),r/(u;)) : w G 0}, 

where co stands for convex hull. In probabilistic terms: S{u, v) is the rate function for joint 
large deviations of (j, r]j). If {u, v) is inside V then we have 

0{u, v) = Su{u, v), t{u, v) = Sv{u, v). 

With slight abuse of notation we shall denote: TTg(^u,v),r{u,v) =■ t^u.v, '^e(u,v),T(u,v) =• 
^e{u,v),T{u,v) =• ^u,v, etc. Clearly, tTu^v can be defined naturally on the boundary of V, in 
that case TTu^i, does not puts zero weight on some of the elements of il. 

We introduce the flux of the conserved quantities. The infinitesimal generator acts on 
the conserved quantities as follows: 

where 

(t){uJi,u]2) ■■= ^ r(wi,u;2;^,W2)(C('^2) - C(^^2)) + C*! 
V'(a;i,u;2) := ^ r(a;i, 0^2; w'^, u;2)('?(^^2) " ??('^2)) + (^2 

(The constants Ci,C2 may be chosen arbitrarily, we will fix them later.) We shall denote the 
expectations of these functions with respect to the canonical measure vr^ ^ by 

$(n, v) := Eu,M, ^{u, v) := E„,„(^/'). (8) 

The following lemma was proved in |12j . 

Lemma 1. Suppose we have a particle system with two conserved quantities and rates satisfying 
conditions Qp and / T^) . Then 

de^ {u{e, t),v{9, r)) = 9,$ {u{e, t),v{9, r)) . 

The first derivative matrix of the fluxes and ^ (with resp. to u, v) will be denoted by 

From Lemma n it follows that D{u,v) is (weakly) hyperbolic, it can be diagonalized in a real 
sense (see \1'2\). We denote the two eigenvalues of D by A and /U, and the corresponding right 
and left eigenvectors by r = (ri,r2)^s = (si,S2)^ and 1 = (/i,/2),m = (mi, 7712): 

Dr = Ar, ID = Al, 
Ds = fis, mD = fim. 
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Although we do not denote it exphcitly, all of these are functions of {u,v). We can assume 

|r| = |s| = 1, 1 • r = 1, m • s = 1. 

The second derivatives of the macroscopic fluxes are denoted by these are symmetric 

two-by-two matrices depending on {u,v). 

2.4 The spectral gap condition 

Let / be a positive integer and {Z,N) integers with Z G [/ min^, Z max N S [/ minr/, / maxry]. 
Expectation with respect to the measure vr^ is denoted by E^_^( • ) . For / : n^z,N ^ ^ let 

i-i 

j = l 

j = l \UJ',LU" J 

L^z N is the infinitesimal generator restricted to the hyperplane ^2'^^, and D^^ is the Dirichlet 
form associated to L'^ ^ (or to its symmetric part). Note, that ^ is defined with free boundary 
conditions. 

We will assume the following additional condition on our models: 

(D) There exists a positive constant W independent of /, Z, N such that for any / : ~^ ^ 

with B^z,Nf = 

B'N,zf<WfD'zMf)- 

Remark. Remark. Presumably (|D|) is true for all (or a large class of) the models satisfying 
conditions (0-(ini)i the techniques of [7j should be suitable to get the desired gap estimates. 
We do not know about any published results covering our case. 

3 Perturbation of the Eulerian hdl 

In ^21 it was proved by the application of Yau's relative entropy method, that under Eulerian 
scaling the local density profiles of the conserved quantities evolve according to the following 
system of partial differential equations: 

dtu + dx^(u,v) = 

(10) 

dtv + dx'^{u,v) = 0. 

This pde is usually a strictly hyperbolic conservation law (i.e. D{u,v) has two distinct real 
eigenvalues) , weak hyperbolicity follows from Lemma Q (see ^2] ) . Since the relative entropy 
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method needs smoothness conditions for the solution of the hmiting equation, the previous 
resuh holds only up to a finite time, till the appearance of the first shock. We also note, 
that jS| gives the first major result about the Eulerian hydrodynamic limit of multi-component 
hyperbolic systems, also wit the application of Yau's method. 

3.1 Formal perturbation 

We will investigate the hydrodynamic behavior of small perturbations of an equilibrium point. 
For that we need to understand the asymptotics of small perturbations of a constant solution of 
(|10j) . One of the perturbation techniques is the so-called method of weakly nonlinear geometric 
optics (see e.g. |S1IS1) which gives the following formal result. 

Fix a point {uo,vo) in V and suppose that this point is strictly hyperbolic, i.e. 

A///, (11) 

at (uqjVq). Suppose {us{t,x),Vs{t,x)) is the solution of the pde (|l()j) with initial conditions 

Me(0,x) = uo + eu*{x), 
VeiO,x) = vo + ev*{x), 

where u*{x),v*{x) are fixed T i-^ R smooth functions. Denote 

ao{x) := 1 • {u*{x),v*{x))'f, := fr^ao{y)dy, 

6o{x) := m • {u* {x) , v* {x))^ , cs := fj6o{y)dy, 



(12) 



and 



ai :=1- (rt$"r,rt^"r)t, aa := 1 • (r^ s, rt ^" s)t, 
6i := m • (s"!" (I)"s, si" ^'"s)''', 62 := m • (r"!" s, ^"s)^ 



(13) 



where 1, m, r, s and are the respective vector- and matrix-valued functions taken at 

Then, according to the formal computations of the geometric optics method, 

MmI ) = ( "0 ) +-(-*---^') ( ) resist..-, t) ( II ) (14) 
as e — > 0, where a and 6 are the solutions of the following Cauchy problems: 

{dta{t,x) + dx {ai ■ ^a{t,x)'^ + csa2(j{t,x)) = 0, 
a{0,x) = ao{x), 



(15) 



and 

r dt6it,x) + dx{bi-l6{t,xf +c^b2 5it,x)) = 0, 
1 5{0,x) = 6o{x). 



(16) 
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Remarks 

1. This result means that a small perturbation of a constant solution of (|1U() is governed by 
the solutions of two decoupled equations (at least, by formal computations). If ai and bi are 
nonzero, then these equations are linear transforms of the Burgers equation. Otherwise the 
respective equations become linear transport equations. It is easy to check, that ai ^ 0, bi ^ 
hold exactly when the point {uo,vq) is genuinely nonlinear, i.e. 

VA • r 7^ 0, V/i • s / 

at iuo,vo). 

2. The geometric optics method is based on series expansion, thus it needs smoothness as a 
condition which could only be true up to a finite time in our case. Surprisingly, this formal 
method gives good approximation of the solutions even after the shocks. In [3j the authors prove 
that the equation (|14j) is valid, in the sense that for any t > the Li-norm of the difference of 
the two sides is bounded by Cte^. In fact, this result is valid for the case if we consider the 
pde UlUf) on T (as we do), on M they have even stronger bounds. 

3.2 The main result 

Our main theorem is a similar result on the microscopic level. We will apply Yau's method, 
thus our results will hold in the regime of smooth solutions, only up to a finite time before the 
first appearance of shocks. 

Suppose, that (no,fo) is a point in the physical domain which is strictly hyperbolic, see 
(|1H) . Let u*(x),v*{x) be smooth real functions on T. Define a{t,x),6{t,x) according to (|12j) . 

(|15() and p6|) . and suppose that they are smooth in T x [0,T]. Fix a small positive 
parameter (3, and suppose that a particle system on fi" satisfying conditions (P | - (|D|) has initial 
distribution for which the density profiles of the two conserved quantities are 'close' to the 
functions uq + n~^u*{-), vq + n~^v*{-). I.e. the profiles are a small perturbation of the constant 
(«o,fo) profile. We also assume, that {uo + n~^u*{x), vo + 'n~^v*{x)) G V holds for every x ^T, 
at least for n > uq. Then, uniformly for < t <T, at time n^^^t the respective density profiles 
will be 'close' to the functions uq + n~^u^'^\t, ■),vq + n~^v^"'\t, •), where 

( :} ) - -f*. - - ^"'*) ( ) + ^(*. - - ( ) ■ 

For the precise formulation of the result we need to introduce some additional notations. We 
will denote by fx^ the true distribution of the system at microscopic time 

:= ;U^xp{ni+^tL-} (18) 

We define the time-dependent reference measure z^p as 
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with defined in ()17() . This measure mimics on a microscopic level the macroscopic 

profiles uq + n~^u^"'\t, •), fo + n~f^v^"'\t, •). We also choose an absolute reference measure 

vr":=nM,<' (20) 

jr'GT" 

which is a stationary measure of our Markov process on The point (uq, Vq) is chosen in a 
way that it lies inside the domain V and 

\uo - + \vo - vi;\ < n-'^ . (21) 

If (no, vq) is inside V, then we may choose (uq , ) = (uq, wq)- By choosing (uq , Vq) inside V we 
get that any probability measure on is absolutely continuous with respect to vr". Condition 
(|21|) ensures that vr" is 'close enough' to /i" in entropy sense, uniformly in t. 

Theorem. Let /? G (0, ^) be fixed. Under the stated conditions, if 

Hif,^\uS) = o{n^-^^), (22) 

then 

i/(/iri^r) = o(ni-2/3), (23) 

uniformly for < t < T . 

The following corollary is a simple consequence of the Theorem and the entropy inequality. 

Corollary. Assume the conditions of Theorem \3.'A Let g : T ^ M 6e a test function. Then for 
any t G [0, T] 



9(-) (0(?^^^^i) -uo) - f g{x) {a{t,x- \nH)ri + 5{t,x- nn^t)si^ dx 
{mi^^^'^i) -vo) - j 9{x) [a{t,x- \nH)r2 + 5{t,x - iJLnH)s^ dx 



n 



^0, 



^0. 



Remarks. 

1. The Theorem states that if the initial distribution of the system is 'close' to ™- relative 
entropy sense then at time n^'^^t it will be close to f". The fact, that 'close' should mean 
o(?i^~^^) can be easily explained, see e.g. jTJ or [T^ . 

2. If instead of condition (|D|) we assume a similar uniform bound on the logarithmic-Sobolev 
constant then our Theorem is valid for /? € (0, ^). 
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4 Proof 



We will assume, that 

{uo,vo) = (0,0), ^>„(0,0) = ^-^(0,0) =0. (24) 

It is easy to see, that we can always reduce the general case to get (|24|) . via some suitable linear 
transformations on (Cv)- Also, with the proper choice of the constants in the definition ((7| we 
can set 

$(0,0) = ^-(0,0) = 0. (25) 

Assumptions (j2l]) imply, that 

I?=(o ly l = rt = (l,0), m = st = (0,l), (26) 

and 

u^"'\t,x)=a{t,x-\n'^t), v^'^\t,x) = 5{t,x- ^nH). (27) 
We introduce the notations 



Clearly, these definitions agree with the definition (|13)) . 
We define the functions a(t,xi,X2),S{t,xi,X2) as 



a{t, xi,X2) ■■= [a2a{t,xi)6{t,X2) + a2(Txit,xi) {6{t, z) - cs)dz + 6{t,X2) 

^ — \ Jo ^ 



2 

(29) 



S{t, xi,X2) := r {b2cr{t,xi)6{t,X2) + b2Sx{t,X2) {a{t, z) - Ca)dz + — a{t, xif 

- A V Jo 2 

The defining partial differential equations ()15p. (|16() of the functions a, 5 are conservation laws, 
thus for any < t < T: 

/ a{z,t)dz = Ca, / S{z,t)dz = cs- 
Jt Jt 

From that it follows that a, 6 are well-defined smooth functions on [0, T] x T x T (i.e. periodic 
in xi and X2) with bounded derivatives. 

4.1 Changing the time-dependent reference measure 

The usual way to prove a result like Theorem l3.2l is to get a Gronwall-type estimate on H{p^\vi): 



Jo 
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via bounding the derivative We will use a slightly different approach, by proving 

a similar estimate for 



Jo 



Here 



n 



vr 



(30) 
(31) 



(32) 



and u^"'\v^'^^ are smooth functions defined as 

5(")(x,t) := u^"\x,t) +n~'^a{t,x - Xn'^t,x - unf^t) 

= a{t, X — Xn^t) + n~^a{t, x — Xn^t, x — fxn^t), 
iy(")(x,t) := v^''\x,t) + n-'^6{t,x- Xn^t,x- unf^t), 

= 6{t,x - pLu^t) + ^(t, X - Xnf^t, x - ^xn'^t) . 
Because of Lemma|^and condition (|22() we have H{^'q\uq) = o(n^~^^) and therefore from H3U() 

H{,,^\1^^) = o(ni-2/5) 
will follow uniformly for < t < T. Using Lemma |21 again we get Theorem 13.21 
Lemma 2. Let /i",i/^,z/" be the measures defined as before, with t G [0,r]. Then 

Hifx^lu^) = o{n^-'^) ^ H{i,^K) = o(ni-2/3). 

Proof. We start with 

H{^,^\v^) - H{f,^\V^) = - log 

By subsections 12.21 and 12.31 we can calculate that 
duj" 

+ (r(n-V"\n-^w(")) - T(n-'^S("), rj-^w^"))) rjj 

-G (^(n-^n("\ n-^^;(")), r(?i-^u("\ n-^r;("))) 

+G (6'(n-^n("),n-^iy(")),r(n-^u("\n~^?y(")))} , 
where, for typographical reasons, we omitted the arguments (t, -) from the functions 



(33) 



From the previous expression via power-series expansion: 



dv^ 



log^(u;) < 0{n^-^^) + Cn-^^ 



jgT" 



G--^^^"^(i,-) 



+ 



+ 



n 



W(i,^; 

n 
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with uniform error terms. Using this with (|33() and the entropy inequahty with respect to 
and t'" the lemma foUows. □ 

We also note, that applying the same arguments as in the proof of Lemma |21 we get that 
from the condition (|22|) 

follows. Since vr" is a stationary measure, 

i7(Mr|vr")<i7(M^|vr") = 0(ni-2'3) (34) 

for all t > 0. 

The proof of the following lemma is a simple application of the entropy inequality with the 
entropy bound jHH). Mind that because of JSU and (f25|l 

Lemma 3. Suppose bi,b2, ■ ■ ■ are real numbers with \bj \ < 1 and stands for either of r]j, 
ipj or . Then 



I ^ b^^^ < Cn 

/On li 



with an absolute constant C. 

In the rest of the paper we prove inequality (jHUJ. 

4.2 Preparatory computations 

We define 



QTl 
^0 



Tn 



nl^e{n-^u^'''^ {t, x) , n-f^v^"^^ {t, x)), 

n^T{n-^u^''\t, x), n-^u(")(i, x)), (35) 



WUln ^UQ,n ^Vq), 



It is easy to check, that the partial derivatives dx0^'^\t,x),dxT^"'\t,x),d^9^'^\t,x),d^T^'^\t,x) 
are uniformly bounded in [0, T] x T. From subsection 12.21 we have 

h 



exp {n-^(^(")(t, ^) - e^)Q + n-P{T(^\t, I) - t^)vj (36) 
-G (n"Pe^^\t, i),n-^T^^\t, i)) + G (n~^ei^\n-^T^''^' ' 



n n 
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Differentiating tlie identity 

and noting that 5t-ff(/i"|7r'") < we get the following bound on 9(if(/i"|z?"): 

n^f-'^tH{^^^\u^) < - [ U^L^ log ^" + n-^+^^dt log dfi'l. 
Integrating with respect to the time: 

Jo Jn^ ^ ' 
We estimate the two terms on the right-hand side separately in the next two subsections. 

4.3 Estimating the first term of (1371) 

From the definitions 



n 



+ (V', - ^{n-Puf\n-Plf;^)) ^) (38) 

+Err5*(t,u;)+Err^(t), 

where 



Err^(t,^) := n'^+^Z^ ^ ('v"0(")(t, ^) - dj^^\t, ^)] 

+ V, (v"?(")(t, ^) - 5,fW(t, ^)) I , (39) 
Err^(t) := n-^+^f ^ U [n~f'uf\n-flfp^) aj(")(t, ^) 

+ {ii-Puf\n-^vf^) d,T^^\t, ^)| . (40) 



We used the (slightly abused) shorthanded notations 

■' n ■' n 

and V" denotes the discrete gradient: 

V'^g{x) ■.= n(^{x+'^)- g{x)^ . 

Using the smoothness of 9^"^^ and r^") and Lemma|21the expectation of the first error term can 
be easily estimated: 

/ |Err?(t,u;)|d^r = 0(n-i+^). (41) 

JO." 
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Using Lemma Hit is easy to see that there exists a smooth function U{u,v) such that 



deU {u{9, r), v{e, r)) = $ {u{e, T),v{e, r)) , drU {u{9, r), v{e, r)) = ^ {u{e, r), v{e, r)) . 
Thus Err2(t) takes the form: 

Err"(i, w) = n^/^-i V d^U ( n'^u^^'Xt, ^) 

from which 

Err2(t) = 0(n-i+2/5), 
uniformly for t G [0,r]. Prom previous bounds we have 

n3/^L"log7;(a;)d/x-ds = 

-1+2/3 r /■ ^ _ ^{n~''uf\n~^vf)) dj^^\s ^ 



(42) 



Jn 



n 







+ 



n 



J J } n 



+0(n-^+^'^) 

In the next step we introduce the block averages. We will denote the block size with / = /(n), 
it will be large microscopically, but small on the macroscopic scale. In the first computations 
we only assume / ^ n^^, the exact order of / will only be determined at the end of the proof, 
after collecting all the error terms. For a local function (j G T") we define its block average 
with 

/-I 



i=0 



n' 







: T M we 


have 


< II^xpIIoo 









Using this with Lemma |21 we can replace ipj in 1)43^ with the the respective block averages: 



n3/^L"log/;(^)d/i^ds 



Jo." 



n 



-1+2/3 



Jn 



jeT" 



+ (^j - ^{n-^uf\n-^vf )) d^T^^\s, ^)| d/^r^^s (44) 
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Finally, using Lemma ^ (the one-block estimate), we replace the block averages (ppi^j by their 
'local equilibrium value': ^{Cpflj) ^^'^ ^(Cji'^j)) respectively: 



/*/ n'f log f 2 {ui)dfiys = 
Jo JU" 




The proof relies on the spectral gap condition (|D|) . It uses the Feynman-Kac formula, 
the Raleigh-Schrodinger perturbation technique and the 'equivalence of ensembles' (see the 
Appendix of [0] for all three). A detailed proof can be found in for the one component case 
which can be easily adapted for our purposes. 

Remark. If instead of the condition (|D|) we have a similar uniform bound on the logarithmic- 
Sobolev constant, then the previous lemma may be strengthened: it holds with the bound 
C{n~^-^fp + l-'). 

4.4 Estimating the second term of (j37|) 

Performing the time-derivation we obtain: 

By the definitions of u^") , v^^'^ and Taylor-expansion we readily get that 

dtu^'^Xt, x) = at{t, X - Xnf^t) - Xn'^a^it, x - Xn'^t) 

— Act^j {t, X — Xn^t, X — fin^t) — fiax2 ^ ~ Xn^t, x — fxn^t) 
+0(n-^). 
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and 

= Xn^axit, X — Xn^t) 
+Xaxi ^ ~ Xn^t, X — fin^t) + Xax2 (t, x — Xn^t, x — fin^t) 
+dx (^o,i<^{t, X — Xn^t)^ + csa2a{t, x — Xn^t)^ 

+dx ^a2cr(t, X — Xn^t){6{t, x — fj,n^t) — cs) + ^a35{t, x — ^n^t)' 

with uniform error terms, i^x^ and 'ax2 are the partial derivatives of xi, X2) with respect 
to the second and third variable.) Adding up these equations and checking the definitions for 
o", (5, a we see that all the significant terms on the right hand side cancel to give: 

+ dx (n2^$(n-^S("), n-^^y("))) = ©(n"^). (47) 

Similarly, 

dtv^""^ + dx (n2^*(n-^5("),n-^iy("))) = ©(n"^). (48) 

From (gZI) and l|l8|) : 

= -n^^Oudx ($(?i-'^5("\n-^^^("))) - n^f^O^dx (^(n-^5("\ n-'^I^^"))) + ©(n"^) 

= -n^ieu^u + e,^u)dxu^''^ - n^ieu^v + e,,^,)dxi^^'^ + o(n-^) 

= -n^ieu^u + r„^ja,S(") - nP{e,^u + r,*„)a,??(") + ©(n"^) 

= -n^$„5^0(") - + 0{n-^) (49) 



In the fourth line we used = 6^ and Lemma ^ To simplify notations, we omitted the 
arguments {t,x) from the functions 9^"'\t^''^\u^'^\v^"'\ and the arguments {n~^u^"'\n~^v^"^) 
from all the partial derivatives of 9, r, ^ with respect to u, v. Similarly, 

= -n'5$„a^^(") - n^^^dxT^''^ + ©(n-'^). (50) 
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Hence from ()46|): 



-n 



+ 



where 



(O - n-^uf^) (n-^Sf \n-^?ff ) ^) (51) 

+Err3(i,ct;), 

Err3(t,^) = i 5] (C, - n-^uf^) b,{t) + (r?, - n'^^jf c,(t) (52) 

jGT" 

and bj{t) and Cj(t) are uniformly bounded constants. Using Lemma |21 we get that 

/ \Eii3it,ui)\dfi'l = Oin~^). (53) 

We can exchange Cj^'Hj with their block-averages Cjj'Hj {^^ the previous subsection) which 
(after the time-integration) gives the following estimate: 

Jo JH" 

+ (Cj - n-Z^nf )) ^„ (n-/3s("),n-/5lf;.")) a,.fW(s, ^) (54) 




+0(n-^ Vn^-^/), 

4.5 Block replacement 

For a function T{u,v) we denote 

Rt('Ui,wi;'U2,W2) := T(ui,fi) - T(n2,f2) - '^uiu2,V2)iui - U2) - T^{U2,V2){VI - V2). 
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Collecting the estimates of the previous subsections we have 

H{^^^\u^)-H{^,^\^;^) 



Jo. 



n 



n 



(55) 



The second derivatives of $ and ^ are bounded thus 

<c({d- n-^u^^\s, ^-)f + {li - n-^v^^\s, ^-)f 



which means that it is sufficient to estimate 

2 



n 



2/3- 



j6T" 



uniformly in t. We estimate the first expression, the other will follow the same way. We denote 



i=0 



g:=g(t,^) = |5:(o-n-/'5(")(t,^) 



bmce x) is uniformly bounded for {t, x) G [0, T] x T, we have 

uniformly in j G T", i € [0, T] and it is enough to estimate 



n 



2/3-1 



Applying the entropy inequality with respect to the time-dependent reference measure and 
using Holder's inequality: 



2/3- 



/„ E (3)'^^" ^ -n2/3-ii7(;,['|i7;^) + i/n^/^-i logE.nexp [^liC^f) , (56) 

for any 7 > 0. P is compact, C is bounded thus there exists a positive constant C such that 

log Bu,v exp ((C - n)y) < Cy^ 

for all (u, v) (z T> and y S M. Thus as a consequence of Lemma El there exists a small 7 > for 
which 



i J] log E,. exp (7/(3)^) 



< 1. 
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Substituting into 



2/3-1 




1 — 



Collecting all the estimates, from (|55p we get 



■t 



n' 



2/3-1 



Choosing / with 




the error term becomes o(l) and the Theorem follows via the Gronwall inequality. (If we have 
the logarithmic-Sobolev condition, and thus a stronger version of LemmalU then / can be chosen 
with <C I <^ to make all the error terms o(l).) 

The proof of Lemma [5] can be found in jllj or |13j . 

Lemma 5. Suppose ^^1,^2) ■ • ■ independent random variables with E^j = for which 



with a positive constant C independent of i and y. Then there exists a small positive constants 
7 depending only on C such that 
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